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We present results on the production of high transverse momentum 7r° and rj mesons in Tv~p 
and 7r~Be interactions at 515 GeV/c. The data span the kinematic ranges 1 < pt < 11 GeV/c in 
transverse momentum and —0.75 < j/cm < 0.75 in rapidity. The inclusive tt" cross sections are com- 
pared with next-to-leading order QCD calculations and to expectations based on a phenomenological 
parton-fcT model. 



PACS numbers: 12.38.Qk, 13.85.Ni 
I. INTRODUCTION 

The study of inclusive single-hadron production at 
large transverse momentum (pr) has been a useful probe 
in the development of perturbative QGD (PQCD) [1, 2]. 
Early in the evolution of the parton model, a departure 
from an exponential dependence of particle production at 
low pt was interpreted in terms of the onset of interac- 
tions between pointlike constituents (partons) contained 
in hadrons. Large pt is a regime where perturbative 
methods have been applied to QCD to provide quantita- 
tive comparisons with data. Such comparisons yield in- 
formation on the validity of the PQCD description, and 
on the parton distribution functions of hadrons and the 
fragmentation functions of partons. 

This paper reports high-precision measurements of 
the production of tt*^ and rj mesons with large px by 
515 GeV/c 7r~ beams. The data were accumulated dur- 
ing the 1990 and 1991-92 fixed-target runs at Fermilab. 
The 7r° production cross sections are compared with next- 
to-leading order (NLO) PQCD calculations [3]. As illus- 
trated in a previous publication [4], our data for both 
inclusive 7r° and direct-photon production are not de- 
scribed satisfactorily by the available NLO PQCD cal- 
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culations, using standard choices of parameters. Sim- 
ilar discrepancies have been observed [5] between con- 
ventional PQCD calculations and other measurements of 
high-pT TJ"" and direct-photon cross sections (see also [6- 
8] ) . The origin of these discrepancies can be attributed 
to the effects of initial-state soft-gluon radiation. Such 
radiation generates transverse components of initial-state 
parton momenta, referred to below as kx [9] . Evidence of 
significant kx in various processes, and a phenomenolog- 
ical model for incorporating its effect on the calculated 
high-pT cross sections, have been extensively discussed 
in Refs. [4, 5]. Recent studies of the photoproduction 
of direct photons at HERA may provide additional in- 
sights [10-14]. In this paper, we follow the phenomeno- 
logical prescription of Ref. [5] when comparing calcula- 
tions with our tt'^ data. 



II. EXPERIMENTAL APPARATUS 

A. Spectrometer 

Fermilab E706 was a fixed-target experiment designed 
to measure the production of direct photons, neutral 
mesons, and associated particles at high-pT [4, 15-17]. 
The spectrometer, designed and built for this experiment, 
was located in the Meson West experimental hall and in- 
cluded a precision charged particle tracking system and a 
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FIG. 1: Plan view of the Meson West spectrometer. The Z- 
axis of the coordinate system for the experiment was oriented 
along the nominal beam axis. 



large acceptance liquid argon calorimeter. Figure 1 shows 
the key elements of the Meson West spectrometer [18]. 

The tracking system consisted of silicon microstrip de- 
tectors (SSDs) in the target region, and multiwire propor- 
tional chambers (PWCs) and straw tube drift chambers 
(STDCs) downstream of a large aperture analysis mag- 
net [17]. The SSD system contained sixteen planes of 
silicon wafers, arranged in eight modules. Each module 
contained two SSD planes, one providing X-view infor- 
mation and the other providing y-view information. Six 
3x3 cm^ SSD planes were located upstream of the tar- 
get and used to reconstruct beam tracks. Two hybrid 
5x5 cm^ SSD planes (25 /im pitch strips in the cen- 
tral 1 cm; 50 /im beyond) were located downstream of 
the target. These were followed by eight 5x5 cm'^ SSD 
planes of 50 //m pitch. The analysis dipole imparted 
a transverse momentum impulse in the horizontal plane 
of « 450 MeV/c to charged particles. Downstream track 
segments were measured by means of four stations of four 
views (XYUV) of 2.54 mm pitch PWCs and two stations 
of eight {AX4Y) layers of STDCs with tube diameters 
1.04 cm (upstream station) and 1.63 cm (downstream 
station). 

In the 1990 run, the target consisted of two 0.8 mm 
thick copper foils followed by two pieces of beryllium. 
The upstream Be piece was 3.7 cm long, while the down- 
stream Be piece was 1.1 cm long. For the 1991-92 run, the 
target was reconfigured to include a liquid hydrogen tar- 
get [19] contained in a 15.3 cm long mylar flask and sup- 
ported in an evacuated volume with beryllium windows 
at each end (2.5 mm thickness upstream and 2.8 mm 
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FIG. 2: The longitudinal distribution of reconstructed ver- 
tices for the 1990 (top) and 1991-92 (bottom) target configu- 
rations. 



thickness downstream). The liquid hydrogen target was 
flanked by two 0.8 mm thick copper disks upstream, and 
a 2.54 cm long beryllium cylinder downstream. Figure 2 
shows the reconstructed vertex position along Z as deter- 
mined by the tracking system for representative samples 
from the 1990 (top) and 1991-92 (bottom) runs. The in- 
dividual target elements are clearly resolved, as are sev- 
eral of the SSDs. 

Photons were detected in a large, lead and liquid-argon 
sampling electromagnetic calorimeter (EMLAC), located 
9 m downstream of the target [20]. The EMLAC had a 
cylindrical geometry with an inner radius of 20 cm and an 
outer radius of 160 cm. It was divided into four mechan- 
ically independent quadrants, which were further subdi- 
vided electronically to create octants. The calorimeter 
had 33 longitudinal cells, read out in two sections: an 11 
cell front section (w 8.5 radiation lengths) and a 22 cell 
back section (w 18 radiation lengths). The longitudinal 
cells consisted of 2 mm thick lead cathodes (the first cath- 
ode was constructed of aluminum) , double-sided copper- 
clad G-10 radial (R) anode boards followed by 2 mm 
thick lead cathodes, and double-sided copper-clad G-10 
azimuthal ($) anode boards. There were 2.5 mm argon 
gaps between each of these layers in a cell. The copper- 
cladding on the anode boards was cut to form strips. 
Signals from corresponding strips on all R (or $) anode 
boards in the front (or back) section were jumpered to- 
gether. The copper-cladding on the radial anode boards 
was cut into concentric strips centered on the nominal 
beam axis. The width of the strips on the first R board 
was 5.5 mm. The width of the R strips on the following 
R boards increased slightly so that the radial geometry 
was projective relative to the target. The azimuthal read- 
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out was subdivided at a radius of 40 cm into inner and 
outer segments, with each inner $ strip subtending an 
azimuthal angle of 7r/192 radians, and outer $ strips cov- 
ering 7r/384 radians. Subdivision of the azimuthal strips 
in the outer portion of the detector improved both the 
position and energy resolution for showers reconstructed 
in this region. It also reduced if!-$ correlation ambi- 
guities from midtiple showers in the same octant of the 
calorimeter. 

The spectrometer also included two other calorimeters: 
a hadronic calorimeter located downstream of the EM- 
LAC within the same cryostat, and a steel and scintilla- 
tor calorimeter positioned further downstream to provide 
instrumentation in the very forward region. The E672 
muon spectrometer, consisting of a toroidal magnet, scin- 
tillators, and proportional wire chambers, was deployed 
immediately downstream of the calorimeters [21-24]. 

The beamline was instrumented with a differential 
Cherenkov counter [25, 26] to identify incident pions, 
kaons, and protons in secondary beams. This helium- 
filled counter was 43.4 m long and was located ~ 100 m 
upstream of the experiment's target. At the end of the 
beamline was a 4.7 m long stack of steel surrounding the 
beam pipe and shadowing the EMLAC to absorb off-axis 
hadrons. A water tank was placed at the downstream end 
of this hadron shield to absorb low-energy neutrons. Sur- 
rounding the hadron shield and neutron absorber were 
walls of scintillation counters (VW) to identify penetrat- 
ing muons. During the 1990 run, there was one wall at 
the upstream end and two walls at the downstream end 
of the hadron absorber. During the 1991-92 run, an ad- 
ditional wall was added to the upstream end. 



B. Trigger 

The E706 trigger selected events yielding high trans- 
verse momentum showers in the EMLAC. The selection 
process involved several stages: beam and interaction 
definitions and the pretriggcr and liigh-p^ trigger require- 
ments. Beam particles were detected using a hodoscope 
consisting of three planes (arranged in X, Y, and U 
views) of scintillator strips, located w 2 m upstream of 
the target region. A beam signal was generated if the 
beam hodoscope registered hits in time coincidence from 
at least two of the three planes. In addition, a BEAmI 
signal was generated if the hits in at least two of the 
hodoscope planes were considered to be isolated. Scin- 
tillation counters centered on the target with a 0.95 cm 
diameter central hole were used to reject interactions ini- 
tiated by particles in the; beam halo. 

Interactions were detected using two pairs of scintil- 
lation counters mounted on the dipole analysis magnet, 
one pair on the upstream side and one pair downstream. 
Each pair had a central hole that allowed non-interacting 
beam particles to pass through undetected. An interac- 
tion was defined as a coincidence between signals from 
at least two of these four interaction counters. A fil- 



ter was used to reject interactions that occurred within 
60 ns of one another to minimize potential confusion in 
the EMLAC due to out-of-time interactions. 

For those interactions that satisfied both the beamI 
and INTERACTION definitions, the pt deposited in various 
regions of the EMLAC was evaluated by weighting the 
energy signals from the fast outputs of the EMLAC R 
channel amplifiers by « sin^j, where 9i is the polar angle 
that the i*'' strip subtended with respect to the nominal 
beam axis. The pretrigger hi signal required that the 
Pt detected in either the inner or the outer R channels 
of any octant was greater than the pretrigger threshold 
value. This signal was issued only when there was no 
evidence in that octant of substantial noise or significant 
Pt attributable to an earlier interaction, and that there 
was no incident beam halo muon detected by the VW. 

Localized trigger groups were formed for each oc- 
tant by clustering the inputs from the R channels into 
32 groups of 8 channels. Each of the adjacent pairs of 
8 channel groups formed a group of 16 strips. If the pt 
detected in any of these groups of 16 was above a speci- 
fied high (or low) threshold, then a LOCAL hi (or LOCAL 
Lo) signal was generated for that octant. If a local hi 
(or LOCAL LO) signal was generated in coincidence with 
the PRETRIGGER HI in the same octant, then a single 
LOCAL HI (or single LOCAL LO) trigger was generated 
for that octant. 

Trigger decisions were also made based upon global en- 
ergy depositions in an octant. A GLOBAL LO signal was 
generated if the total pt in an octant exceeded a low 
global threshold value. The LOCAL(g)GLOBAL LO trig- 
ger required a coincidence of the GLOBAL LO signal with 
PRETRIGGER HI and LOCAL LO signals from the same 
octant. The local LO requirement was included to sup- 
press spurious global triggers due to coherent noise in the 
EMLAC. 

The SINGLE LOCAL LO and LOCAL0GLOBAL LO trig- 
gers were prescaled to keep them from dominating the 
data sample. Prescaled samples of beam, interaction, 
and pretrigger events were also recorded. Further de- 
tails concerning the E706 trigger can be found else- 
where [15, 17, 27, 28]. 



III. ANALYSIS METHODS 

This paper presents results for 7r° and rj production by 
515 GeV/c 7r~ beam on beryllium and liquid hydrogen 
targets. The data sample used in this analysis corre- 
sponds to an integrated luminosity of 7.5 pb~^ of 7r~Be 
data (6.1 pb~^ from the 1990 run plus 1.4 pb~^ from the 
1991-92 run) and 0.23 pb-^ of n'p data. The following 
subsections describe the data analysis procedures and the 
methods used to correct the data for losses due to inefla- 
ciencics and selection biases. Additional details may be 
found in several of our previous papers [15, 17, 20]. 
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A. Reconstruction 

Two major aspects of the event reconstruction proce- 
dure were particle track and calorimeter shower recon- 
struction. The charged-track reconstruction algorithm 
produced track segments upstream of the magnet using 
information from the SSDs, and downstream of the mag- 
net using information from the PWCs and STDCs. These 
track segments were projected to the center of the magnet 
and linked to form the final reconstructed tracks, whose 
calculated charges and momenta were used in physics 
analyses and to determine the location of the primary 
interaction vertex. The charged track reconstruction is 
described in detail in Ref. [17]. 

The readout in each EMLAC quadrant consisted of 
four regions: left R and right R, (radial strips of each oc- 
tant in that quadrant), and inner $ and outer Strip 
energies from clusters in each region were fit to the shape 
of an electromagnetic shower as determined from Monte 
Carlo simulations and isolated-shower data. These fits 
were used to evaluate the positions and energies of the 
peaks in each region. Shower positions and energies were 
obtained by correlating peaks of approximately the same 
energy in the R and $ regions within the same half oc- 
tant (more complex algorithms were used to handle con- 
figurations with overlapping showers in either the R or 
$ regions). The EMLAC readout was also subdivided 
longitudinally into front and back sections. This lon- 
gitudinal segmentation provided discrimination between 
showers generated by electromagnetically or hadronically 
interacting particles. An expanded discussion of the EM- 
LAC reconstruction procedures and performance can be 
found in Ref. [20]. 
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FIG. 3: 77 mass distributions in the region of the tt" for 
several pr bins. Curves are overlayed for those pr bins where 

the background to the signal was determined using a fitting 
procedure rather than sideband subtraction. 
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FIG. 4: 77 mass distributions in the region of the r) for several 
Pt bins. 



B. Event selection and meson signals 

Events contributing to measurements of cross sections 
were required to have reconstructed vertices within the 
fiducial volume of the Be or H2 targets. Both 7r° and 
7] mesons were reconstructed via their 77 decay modes. 
Only those showers which deposited at least 20% of their 
energy in the front part of EMLAC were considered as 
7 candidates, to reduce the background due to showers 
from hadronic interactions. Only those 77 combinations 
with energy asymmetry A^^y = \E-y-^ 1/ (^71 +-^72) < 

0.75 were considered in order to reduce micertainties due 
to low energy photons. Photons were required to be re- 
constructed within the fiducial region of the EMLAC to 
exclude areas with reduced sensitivity. In particular, re- 
gions of the detector near quadrant boundaries (which 
abutted steel support plates), the central beam hole, the 
outer radius of the EMLAC, and octant boundaries were 
excluded. In addition, 77 combinations were restricted to 
the same octant to simplify the trigger analysis. A sim- 
ple ray-tracing Monte Carlo program was employed to 
determine the correction for these fiducial requirements. 

Signals have been corrected for the 25% loss due to the 



energy asymmetry cut and the branching fractions for the 
77 decay modes [29] . The correction for losses due to the 
conversion of one or both of the photons into e+e^ pairs 
was evaluated by projecting each reconstructed photon 
from the event vertex to the reconstructed position in 
the EMLAC. The radiation length of material traversed 
was calculated based upon detailed descriptions of the de- 
tectors encountered. The photon conversion probability 
was evaluated and used to account for conversion losses. 
The average correction for conversion losses was 1.19 for 
the Be target in the 1990 run (1.16 in 1991-92) and 1.24 
for the p target. A full event simulation (described be- 
low) was employed to correct for other effects including 
reconstruction smearing and losses. 

The 77 invariant mass distributions in the 7r° and rj 
mass regions for photon-pairs that satisfied the above re- 
quirements are shown in Figs. 3 and 4 for several 77 pt 
intervals. A tt" candidate was defined as a combination 
of two photons with invariant mass, M^^, in the range 
100 MeV/c^ < M^.^ < 180 MeV/c^. An rj candidate 
was defined as a two-photon combination in the range 
450MeV/c2 < M^^ < eSOMeV/c^. Combinatorial back- 
ground under the tt^ and 77 peak regions was evaluated as 
follows. Sideband regions were defined to cover a mass 
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range equivalent to that in the tt^ and rj peak regions. 
Distributions from these side bands were subtracted from 
the distributions within the tt" and ry mass ranges to ob- 
tain the respective signals. This technique is appropriate 
as long as the combinatorial background depends approx- 
imately linearly on M^j. The combinatorial background 
shape is not linear at low pt (below « 2 GeV/c), and a 
fitting procedure was used to evaluate this background. 
The 77 mass distributions at low px were fit using Gaus- 
sians for signal, and second and third order polynomials 
for the background. The combinatorial background in 
the peak regions was determined from the resultant fit 
parameters, and subsequently subtracted from the total 
entries within the peak. 



C. Trigger Corrections 

Trigger corrections were evaluated on an event-by- 
event basis using the measured efficiencies of the trigger 
groups responsible for the formation of a given trigger. 
For example, the single local hi trigger corrections 
were based upon the efficiencies of the 32 groups of 16 
in the triggering octant. These efficiencies were evalu- 
ated as functions of the pr reconstructed within the trig- 
ger group, using data samples that were unbiased with 
respect to the trigger group. From these efficiencies, a 
trigger probability was defined, P = 1 — n(l where 
Pi is the efficiency of the i*'' trigger group in the oc- 
tant. The inverse of this probability was applied as a 
trigger weight to each meson candidate. Meson candi- 
dates with trigger probabilities of P < 0.1 were excluded 
from further consideration to avoid excessively large trig- 
ger weights. The correction for losses from this require- 
ment was determined from Monte Carlo, and absorbed 
into the reconstruction efficiency. 

The cross sections presented in this paper utilize the 
results from a combination of triggers. Going from low 
to high Pt, the triggers were: INTERACTION, PRETRIG- 

GER, HI, LOCAL(g)GLOBAL LO (1990 run) Or SINGLE LOCAL 

LO (1991-92 run), and single local hi. The transi- 
tion points chosen between the lower and higher thresh- 
old triggers were determined by comparing the fully cor- 
rected results from each trigger, and were different for 
7r° and i] mesons, and also depended on rapidity. For 
Pt> 4.0 GeV/c, the single local hi trigger was used 
exclusively for both tt" and rj mesons. 



D. Beam halo muon rejection 

Spurious triggers were produced by muons in the 
beam halo that deposited energy in the electromagnetic 

calorimeter in random coincidence with an interaction in 
the target. Particularly in the outer regions of the EM- 
LAC, such energy depositions can produce background 
at low 77 mass values due to the occasional splitting 
of the muon-induced showers into two closely-separated 
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FIG. 5: Eflfect of muon rejection requirements on the invari- 
ant mass distribution in the Tr^-mass region, for candidate 77 
pairs with 3 < pt < 3.5 GeV/c (left) and 7 < pr < 10 GeV/c 
(right). 



photon candidates. The pretrigger logic used signals from 
the VW to reject events associated with such muons in 
the beam halo. The off-line analysis employed expanded 
requirements on the latched VW signals, and imposed 
requirements upon the direction of reconstructed show- 
ers, the shower shape, and the total px imbalance in the 
event to futher reduce this background. The effects of 
these off-line rejection requirements on the 77 invariant 
mass distribution are shown in Fig. 5, for 77 pairs with 
3 < Pt < 3.5 GeV/c (left) and 7 < pr < 10 GeV/c 
(right). The rejection requirements completely eliminate 
the large muon-induced background in the high-pT bin 
while having very little effect on the signal. A more de- 
tailed description of these requirements can be found in 
Ref. [15]. 



E. Detector simulation 

The Meson West spectrometer was modeled via a de- 
tailed GEANT [30] simulation. A preprocessor was used 
to convert GEANT information into the hits and strip en- 
ergies measured by the various detectors. The preproces- 
sor simulated hardware effects, such as channel noise and 
gain variations. Monte Carlo generated events were then 
processed through the same reconstruction software used 
to analyze the data. This accounted for inefficiencies and 
biases in the reconstruction algorithms. Reconstruction 
inefficiencies for tt^ and 77 mesons were relatively small 
over most of the kinematic range. More information on 



6 



> 3.5 GeV/c rTi 




Pj > 6.5 GeV/c 

J 




JtBeat515GeV 
-0.75 < y^_,, < 0.75 



0.20 0.22 0.24 

yyMass (GeV/c^) 



s 

o 

S " 4.0 < p.;. < 5.5 GeV/c 



71 Beat515GeV/c 
-0.75 <y,„< 0.75 




1\ 



\L 



0.0 0.25 0.5 0.75 1.0 



0.6 0.8 1.0 

Energy Asymmetry 



FIG. 6: Comparison between data (•) and the Monte Carlo 
(histogram) for 77 combinations in the 7r° mass region for two 
requirements on minimum pr- The distributions have been 
normahzed to unit area. 



the detailed simulation of the Meson West spectrometer 
can be found in Rcf. [15]. Wc employed single particle 
distributions, reconstructed data, and the herwig [31] 
physics generator as inputs to the GEANT simulations. 
HERWIG calculations of tt" and 77 spectra were weighted 
in pT and rapidity using our measured results in an iter- 
ative fashion so that the final corrections were based on 
the data distributions rather than on the output of the 
physics generator. 

Spectral effects were particularly important to the cal- 
ibration of the EMLAC's energy response [20]. The cal- 
ibration of the energy response was based on the recon- 
structed masses of tt*^ mesons in the 77 decay mode. 
The steeply falling 7r° pt spectrum, combined with the 
calorimeter's resolution, results in a small offset (« 1%) 
in the mean reconstructed photon energies. We ac- 
counted for this offset, and for potential biases in the 
calibration procedure, by calibrating the simulated EM- 
LAC in the same manner as the real detector. We also 
employed the simulation to evaluate the mean correction 
(as a function of photon energy) for energy deposited in 
the material upstream of the EMLAC. The impacts of 
detector resolution on the energy scale calibration and 
on the TT*^ and rj production spectra were incorporated in 
the overall reconstruction efficiency corrections. 

Figure 6 compares the 77 mass spectra in the tt^ and 
7] mass regions to the simulation for two different mini- 
mum Pt cutoffs. In addition to providing evidence that 
the Monte Carlo simulated the EMLAC resolution well, 
the agreement between the levels of combinatorial back- 
ground also indicates that the Monte Carlo provided a 



FIG. 7: Comparison of the background subtracted 77 energy 
asymmetry distribution in the 7r° si gnal region in data (•) and 
the Monte Carlo (histogram). Shown arc the comparisons for 
two Pt intervals, 4.0 < pT < 5.5 GeV/c and 5.5 < pr < 
7.0 GeV/c. The distributions have been normalized to unit 
area. 



reasonable simulation of the underlying event structure. 
Figure 7 shows a comparison between the Monte Carlo 
simulation and the data for the background subtracted 
77 energy asymmetry in the tt*^ signal region, for two pt 
intervals. This figure illustrates that the simulation ac- 
curately describes the losses of very low-energy photons. 



F. Normalization 

Electronic scalers that counted signals from the beam 
hodoscope, interaction counters, and beam hole counters 
were used to evaluate the number of beam particles inci- 
dent on the target. Other scalers logged the state of the 
trigger and of components of the data acquisition system. 
Information from these scalers was used to determine the 
number of beam particles that traversed the spectrome- 
ter when it was ready to record data. This number was 
corrected for multiple occupancy in the beam (« 3%), 
the absorption of beam in the target material (w 6% for 
the Be target and « 3% for p), and for the /x~ content of 
the beam which was measured to be « 0.5% [28]. 

The normalization of the low pt 7r° cross section was 
independently verified using events from the prescaled 
beam and interaction trigger samples. In these samples, 
the absolute normalization is obtained simply by count- 
ing the recorded events selected by these triggers. For 
these low pt events, the normalization as determined via 
the scalers and via event counting techniques agreed to 
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FIG. 8: Combined relative systematic uncertainty for 7r° pro- 
duction as a function of pr (solid line) . Also shown are the 
contributions from the various sources of systematic uncer- 
tainty. 



FIG. 9: Ratio of production cross sections by 515 GeV/c tt" 
beam on Be from the 1991-92 fixed target run to those ob- 
tained from the 1990 run. The error bars reflect only statis- 
tical uncertainties for pr > 2.0 GeV/c. The dashed line is a 
fit to the ratio of tt" cross sections. 



3% accuracy. 

An analysis of negative secondary beam production by 
800 GeV/c primary protons indicates a small (« 1%) K~ 
component in the incident 515 GcV/c beam. The meson 
production cross sections were corrected for this com- 
ponent under the assumption that meson production by 
K~ beam is half that of tt" beam [25, 32]. 



G. Summary of systematic uncertainties 

The principal contributions to the systematic uncer- 
tainty arose from the following sources: calibration of 
photon energy response, 7r° and r] reconstruction effi- 
ciency and detector-resolution unsmearing, and overall 
normalization. The relative systematic uncertainty for 
7r° production is shown as a function of in Fig. 8. In- 
cluded in the figure are curves showing the contributions 
from the major sources of systematic uncertainty. Other 
sources of uncertainty which contribute at the 1-2% level 
include: background subtraction, incident beam contam- 
ination, beam halo muon rejection, geometric acceptance, 
photon conversions, trigger response, and vertex finding. 
The total systematic uncertainty is calculated by combin- 
ing in quadrature all the individual uncertainties. The 
uncertainties for r] production are similar to those for 7r° 
production, except for the imcertainty in the trigger re- 
sponse, which is ~ 5% at low px- The actual systematic 
uncertainties are quoted in the appropriate tables of cross 
sections, except at low pT, where, the statistical and sys- 
tematic uncertainties have been combined because of the 



large correlation between them. 

A cross-check of the systematics can be obtained by 
comparing the results for meson production from the 
1991-92 data sample to the results from the 1990 data 
sample for independent analyses. The ratio of the tt" 
and r] production cross sections from the 1991-92 run to 
the 1990 run is shown in Fig. 9. A linear fit to the ratio 
of tt" cross sections yields 0.982 ± 0.006, which is well 
within the systematic uncertainties described earlier in 
this section. 

The secondary pion beams were determined to have a 
mean momentum of 515 ± 2 GeV/c with an estimated 
halfwidth of w 30 GeV/c. This momentum spread in- 
troduces a small uncertainty (~ 5%) in comparisons of 
theory with data. 



IV. RESULTS AND DISCUSSION 

A. 7r° production 

The inclusive 7r° cross sections per nucleon versus pr 
are shown in Fig. 10 for 515 GeV/c tt^ beams incident 
upon beryllium and liquid hydrogen targets. Note that 
the cross section is in units of pb/(GeV/c)^ per nucleon 
for the Be target, and in nb/(GeV/c)^ for the p target. 
The measurements on beryllium were obtained using the 
combined statistics from the 1990 and 1991-92 runs. Be- 
cause of the steeply falling spectra, the data are plotted 
at abscissa values that correspond to the average values 
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FIG. 10: Invariant difFcrcutial cross sections (per nucleon) for 
tt" production as a function of px in Tr~p and ir~ Be interac- 
tions at 515 GeV/c. Cross sections have been averaged over 
the rapidity range, —0.75 < ycm < 0.75. The error bars rep- 
resent the combined statistical and systematic uncertainties. 



FIG. 11: Invariant differential cross sections (per nucleon) for 
tt'^ production as a function of pr in ■K~p and n~ Be interac- 
tions at 515 GeV/c. Cross sections have been averaged over 
the rapidity range, —0.75 < j/cm < 0.75. The error bars have 
statistical and systematic uncertainties added in quadrature. 
Overlayed on the data are NLO PQCD results. 



of the cross section in each pT bin, assuming local expo- 
nential pt dependence [33]. The cross sections are also 
tabulated in Tables I through IV. 

In Fig. 11, the measured inclusive 7r° cross sections 
arc compared to NLO PQCD results [3] using GRV [34] 
parton distributions, KKP [35] fragmentation functions, 
and factorization scales of n = Pt/^, Pt, and 2pT (the 
renormaUzation and fragmentation scales have been set 
equal to the value of the factorization scale). The PQCD 
calculations for the Be target were adjusted to account 
for nuclear effects using the HIJING Monte Carlo calcula- 
tion [36]. The large scale dependence in the calculations 
is insufficient to raise the predicted cross sections up to 
the values indicated by the data. 

These discrepancies have been interpreted [4, 5] as aris- 
ing from additional soft-gluon emission in the initial state 
that is not included in the NLO calculation, and which re- 
sults in sizeable parton kr before the hard collision (for a 
different perspective, sec the discussion in Ref. [8]). Soft- 
gluon (or kx) effects arc expected in all hard-scattering 
processes, such as the inclusive production of jets, high- 
er mesons, and direct photons [37-40]. The CoUins- 
Soper-Sterman resummation formalism [41] provides a 
rigorous basis for understanding these radiative effects, 
and there have been several recent efforts to derive resum- 
mation descriptions for the inclusive direct-photon [42- 
46], jet [47], and dijet cross sections [48 50]. The calcu- 
lation of Ref. [42] for inclusive direct-photon production, 
which includes the effects of soft-gluon resummation near 
the kinematic threshold limit {xt = 2pt/\/s — > I), has 
a far smaller sensitivity to scale, compared to NLO cal- 



culations, and provides cross sections close to those of 
NLO calculations with a scale of /x = pt/2. Also, for 
our energies, the calculations of Ref. [45, 46], which si- 
multaneously treat threshold and recoil effects in direct- 
photon production, yield a substantially larger cross sec- 
tion than the NLO result. However, no such calcula- 
tions are available for inclusive meson production. In 
their absence, we use a PQCD-based model that incor- 
porates transverse kinematics of initial-state partons to 
study the consequences of additional kx for high-pr pro- 
duction processes. 

Because the unmodified PQCD cross sections fall 
rapidly with increasing pt, the net effect of the "fcr 
smearing" is to increase the expected yield. Modified 
parton kinematics have been implemented in a Monte 
Carlo calculation of the leading-order (LO) cross sections 
for high-pT particle production [51], with the kr distri- 
bution for each of the incoming partons represented by a 
Gaussian with one adjustable parameter ((fcr))- Unfor- 
tunately, no such program is available for NLO calcula- 
tions, and so we approximate the effect of kr smearing 
by multiplying the NLO cross sections by the correspond- 
ing LO fcT-cnhanccmcnt factors. Admittedly, this proce- 
dure involves a risk of double-counting since some of the 
fcT-enhancement may already be contained in the NLO 
calculation. However, we expect such double-counting 
effects to be small. 

In the calculation of the LO fcT-enhancement factors 
we employ (kx) values representative of those found from 
comparisons of kinematic distributions in data involv- 
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FIG. 12: Invariant difFcrcutial cross sections (per nucleon) for 
tt" production as a function of px in Tr~p and 7r~ Be interac- 
tions at 515 GeV/c. Cross sections have been averaged over 
the rapidity range, —0.75 < i/cm < 0.75. The error bars have 
statistical and systematic uncertainties added in quadrature. 
The curves represent the /cr-enhanced NLO QCD calculations 
for {fcT)=1.4 GeV/c. 



FIG. 13: Invariant cross sections per nucleon for iv" produc- 
tion in Tr~Be interactions at 515 GeV/c. Cross sections are 
shown versus j/cm for several intervals in pr- The curves 
represent fcT-cnhanced NLO QCD caleulations for {A:t)=1.4 
GeV/c. The error bars have statistical and systematic uncer- 
tainties added in quadrature. 



ing production of high-mass 77, 771", and 7r°7r° systems 
with calculations relying on the same LO program (see 
Rcfs. [4, 5] for further details). For these comparisons, 
we used the LO versions of the GRV parton distributions 
and (where appropriate) KKP fragmentation functions, 
and an average transverse momentum of 0.6 GeV/c for 
the 7r° mesons relative to the fragmenting parton direc- 
tion (varying this parameter in the range 0.3—0.7 GeV/c 
does not affect our conclusions) [52, 53]. 

Comparisons of the fcr-enhanced calculations with 
data at 515 GeV/c arc displayed in Fig. 12, indicating 
good agreement for the chosen {kx) value. Figure 13 
shows the ■K~Be cross sections at 515 GeV/c versus ra- 
pidity, for several intervals in pT- The shapes and nor- 
malizations of the data are in good agreement with the 
fcr-enhanced calculations. 

In a previous publication [15], this experiment reported 
results for tt" production in pp and pQc interactions at 
800 and 530 GeV/c. This offers the opportunity to 
compare production cross sections by incident 7r~ and p 
beams. In Fig. 14, the ratio of invariant cross sections for 
7r° production by 515 GeV/c 7r~ and 800 GeV/c p beam 
is shown as a function of xt (comparing results at ap- 
proximately the same incident momentum per valence 
quark). Both theoretical and experimental uncertainties 
are reduced in the ratio allowing, in principle, a more 
sensitive test of the calculations. In the figure, the data 
are compared to conventional {{kT)=0) and /cT-cnhanced 
NLO results using KKP fragmentation functions. The 



D 20 



1 1 1 1 1 1 1 1 1 1 
7c"Be->7t''Xat515GeV/c 






pBe 7c"x at 800 GeV/c 








NLOTheory (n = p^/2) 








GRV pdf(ji" Beam) 








CTEQ4Mpdf(pBeam) 








KKP fragmentation 








fl.40GeV/c(515GeV/c) 

---- <u= ,-- 1 

^ [1.10 GeV/c (800 GeV/c) i 








<kj> = 0.0 GeV/c \ 




, ♦ 

\ -0.75 <y^^< 0.75, 515 GeV/crt beam 


-1 .0 < y^,^ < 0.5, 800 GeV/c p beam 


Stat uncertainties only 

1 .... 1 .... 1 .. . 



I ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' 1 

0.1 0.2 0.3 0.4 0.5 



FIG. 14: Ratio of 515 GeV/c tt" to 800 GeV/c p beam 7r° 
production cross sections as a function of xt, compared to 
conventional and fcr-enhanced NLO QCD results. 
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FIG. 15: Ratios of 7] to tt" invariant cross sections at 515 
GeV/c, as functions of pr (top), and of j/cm (bottom) for two 
Pt ranges. The error bars reflect only the statistical uncer- 
tainties. 



fcr-enhanced theory accommodates the data better than 
the conventional theory. 



B. r? production 

Cross sections for inclusive t] production are tabulated 
in Tables V through VIII. Theoretical descriptions of 
?7-meson production differ from the tt^ case primarily be- 
cause of differences in the fragmentation of partons into 
the particles of interest. To investigate this effect, we 
present rj/ir'^ relative production rates as functions ofpr 
and 2/cm (for two px ranges) in Fig. 15. We see no sig- 
nificant dependence in this ratio — the data average to 
a value of 0.48 ± 0.01 (statistical). 



V. SUMMARY 

The invariant cross sections for tt*^ and ry production 
have been measured for Tr~p and Tr~Be collisions at 



515 GeV/c as functions oipx and j/cm, over the kinematic 
range 1< pt < H GeV/c, and -0.75 < ycm < 0.75. The 
TT^ cross sections agree with fcr-cnhanced NLO QCD cal- 
culations. The measured tj/tt^ production ratio, which 
provides information about the relative fragmentation of 
partons into these mesons, is 0.48 ± 0.01 with little de- 
pendence on Pt or rapidity. 
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APPENDIX: TABULATED CROSS SECTIONS 



TABLE I: Invariant differential cross section (^Eda/d^p) per 
nuclcon for the inclusive reaction 7r~Be tt^X at 515 GeV/c, 
averaged over the rapidity interval —0.75 < j/cm < 0.75. 



TABLE II: Invariant differential cross section (Eda / d^p) for 

the inclusive reaction -K^p tt'^X at 515 GeV/c, averaged 
over the rapidity interval —0.75 < j/cm < 0.75. 



Pt 
(GeV/c) 


TT-p @ 515 GeV/c 
(/.b/(GeV/c)^) 


1.00 - 1.50 


161 ± 57 


1.50 - 2.50 


5.2 ± 4.1 


(nb/(GeV/c)^) 


2.50 - 2.75 


455 ± 34 ± 52 


2.75 - 3.00 


170 ± 14 ± 19 


3.00 - 3.25 


96.3 ± 8.3 ± 10.9 


3.25 - 3.50 


34.3 ± 3.7 ± 3.8 


3.50 - 3.75 


18.5 ± 2.6 ± 2.1 


3.75 - 4.00 


9.8 ± 1.4 ± 1.1 


4.00 - 4.25 


4.67 ± 0.12 ±0.51 


4.25 - 4.50 


2.586 ± 0.085 ± 0.28 


4.50 - 4.75 


1.303 ± 0.048 ± 0.14 


(pb/(GoV/c)-) 


4.75 - 5.00 


715 ± 33 ± 79 


5.00 - 5.25 


386 ± 24 ± 43 


5.25 - 5.50 


213 ± 16 ± 24 


5.50 - 5.75 


137 ± 13 ±15 


5.75 - 6.00 


69.7 ± 9.3 ±7.8 


6.00 - 6.50 


31.8 ± 4.1 ±3.6 


6.50 - 7.00 


11.2 ± 2.3 ± 1.3 


7.00 - 8.00 


3.69 ± 0.86 ± 0.43 


8.00 - 10.00 


0.32 ± 0.19 ±0.04 



Pt 


7r Be 


@ 515 GeV/c 
/(GeV/c) = ) 


(GeV/c) 


(pb 


1.00 - 


1.20 


255 ± 


35 


1.20 - 


1.40 


97 ± 


14 


1.40 - 


1.60 


39.5 ± 


5.9 


1.60 - 


1.80 


19.8 ± 


2.5 


1.80 - 


2.00 


9.5 ± 


1.2 


2.00 - 


2.20 


3.64 ± 


0.36 ± 0.43 


2.20 - 


2.30 


2.112 ± 


0.053 ±0.25 


2.30 - 


2.40 


1.460 ± 


0.038 ±0.17 


(nb/(GeV/c)^) 


2.40 - 


2.50 


955 ± 


27 ± 111 


2.50 - 


2.60 


673 ± 


17 ± 78 


2.60 - 


2.70 


496 ± 


14 ± 57 


2.70 - 


2.80 


315.4 ± 


9.6 ± 36.1 


2.80 - 


2.90 


243.8 ± 


7.2 ± 27.8 


2.90 - 


3.00 


170.2 ± 


4.3 ± 19.3 


3.00 - 


3.10 


127.5 ± 


3.3 ± 14.4 


3.10 - 


3.20 


90.7 ± 


2.6 ± 10.2 


3.20 - 


3.30 


72.1 ± 


1.9 ±8.1 


3.30 - 


3.40 


48.3 ± 


1.3 ± 5.4 


3.40 - 


3.50 


35.38 ± 


1.00 ± 3.9 


3.50 - 


3.60 


26.53 ± 


0.65 ± 2.9 


3.60 - 


3.70 


19.73 ± 


0.50 ± 2.2 


3.70 - 


3.80 


14.88 ± 


0.36 ± 1.6 


3.80 - 


3.90 


11.15 ± 


0.30 ± 1.2 


3.90 - 


4.00 


8.28 ± 


0.11 ± 0.91 


4.00 - 


4.10 


6.220 ± 


0.039 ± 0.69 


4.10 - 


4.20 


4.770 ± 


0.034 ± 0.53 


4.20 - 


4.30 


3.604 ± 


0.027 ± 0.40 


4.30 - 


4.40 


2.825 ± 


0.022 ± 0.31 


4.40 - 


4.50 


2.166 ± 


0.019 ± 0.24 


4.50 - 


4.60 


1.654 ± 


0.015 ± 0.18 


4.60 - 


4.70 


1.277 ± 


0.013 ± 0.14 


(pb/(GoV/c)-) 


4.70 - 


4.80 


990 ± 


11 ± 108 


4.80 - 


4.90 


756.4 ± 


9.5 ± 83.0 


4.90 - 


5.00 


605.8 ± 


8.5 ±66.6 


5.00 - 


5.10 


475.9 ± 


7.2 ± 52.4 


5.10 - 


5.20 


366.0 ± 


6.0 ±40.3 


5.20 - 


5.30 


298.1 ± 


5.4 ± 32.9 


5.30 - 


5.40 


237.8 ± 


4.8 ± 26.3 


5.40 - 


5.50 


191.4 ± 


4.3 ± 21.2 


5.50 - 


5.60 


151.3 ± 


3.7 ± 16.8 


5.60 - 


5.70 


116.5 ± 


3.2 ± 13.0 


5.70 - 


5.80 


91.0 ± 


2.8 ± 10.2 


5.80 - 


5.90 


77.2 ± 


2.6 ±8.6 


5.90 - 


6.00 


57.3 ± 


2.2 ±6.4 


6.00 - 


6.25 


41.3 ± 


1.2 ± 4.7 


6.25 - 


6.50 


25.27 ± 


0.90 ± 2.9 


6.50 - 


6.75 


15.27 ± 


0.67 ± 1.7 


6.75 - 


7.00 


8.71 ± 


0.50 ± 1.0 


7.00 - 


7.50 


4.48 ± 


0.25 ± 0.52 


7.50 - 


8.00 


1.30 ± 


0.13 ± 0.16 


8.00 - 


9.00 


0.366 ± 


0.049 ± 0.045 


9.00 - 


10.00 


0.044 ± 


0.022 ± 0.006 


10.00 - 


12.00 


0.0031 ± 


0.0031 ± 0.0004 
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TABLE III: The averaged invariant differential cross section (Eda/d^p) per nucleon as a function of rapidity and pr for the 
inclusive reaction 7r~Be — » tt^X at 515 GeV/c. 





Vera 




Pt 


(GeV/c) 














1.00 - 1.50 


1.50 - 2.00 


2.00 - 2.50 




2.50 - 3.00 










(Mb/(GeV/c)2) 


(/:ib/(GeV/c) = ) 


(/*b/(GeV/c) = ) 




(nb/(GcV/c)2: 


) 


-0, 
-0, 


.750 
.625 


- -0.625 
0.500 


132 ± 30 


13.3 ± 3.7 


1.84 ± 0.86 ± 0.22 




373 ± 26 ± 43 
400 ± 24 ± 46 




-0, 


.500 


0.375 


175 ± 29 


18.4 ± 4.4 


2.63 ± 0.13 ± 0.31 




366 ± 24 ± 42 




— 0. 


.375 


- -0.250 


2.48 ± 0.10 ± 0.29 




336 ± 16 ± 38 




— 0. 


.250 


- -0.125 


148 ± 26 


23.2 ± 4.0 


2.538 ± 0.084 ± 0.30 




375 ± 16 ± 43 




— 0, 


.125 


- 0.000 


2.522 ± 0.068 ± 0.30 




404 ± 15 ± 46 







,000 


- 0.125 


171 ± 31 


19.6 ± 4.0 


2.506 ± 0.064 ± 0.29 




409 ± 13 ± 47 




0. 


,125 


- 0.250 


2.551 ± 0.063 ± 0.30 




392 ± 12 ± 45 




0. 


,250 


- 0.375 


124 ± 22 


24.6 ± 4.4 


2.451 ± 0.067 ± 0.29 




393 ± 12 ± 45 




0. 


,375 


- 0.500 


2.360 ± 0.070 ± 0.28 




354 ± 11 ±41 




0, 


,500 


- 0.625 


144 ± 24 


14.8 ±3.2 


2.407 ± 0.074 ± 0.28 




368 ± 12 ±42 




0, 


,625 


- 0.750 


2.083 ± 0.071 ± 0.24 




334 ± 12 ± 38 










3.00 - 3.50 


3.50 - 4.00 


4.00 - 4.50 




4.50 - 5.00 










(nb/(GeV/c)^) 


(nb/(GeV/c)^) 


(nb/(GeV/c)^) 




(nb/(GeV/c)^) 


-0, 


,750 


- -0.625 


60.8 ± 5.7 ±6.8 


10.38 ± 0.97 ± 1.1 


2.198 ± 0.071 ±0.24 


0, 


.614 ± 0.022 ± 0, 


,067 


-0, 


,625 


- -0.500 


61.0 ± 5.8 ± 6.8 


11.59± 0.95 ± 1.3 


2.766 ± 0.057 ±0.30 


0, 


.770 ± 0.019 ± 0, 


,084 


-0 


,500 


- -0.375 


73.3 ± 4.1 ± 8.2 


14.27± 0.46 ± 1.6 


3.309 ± 0.044 ±0.36 


0, 


.863 ± 0.016 ±0, 


,095 


— 


,375 


0.250 


68.6 ± 3.3 ± 7.7 


14.96 ± 0.57 ±1.7 


3.668 ± 0.041 ± 0.40 





,955 ± 0.017 ± 0, 


,10 


-0, 


,250 


0.125 


77.1 ± 2.9 ±8.7 


17.31 ± 0.65 ± 1.9 


4.067 ± 0.038 ± 0.45 


1 


,117 ± 0.018 ± 0, 


,12 


-0, 


,125 


- 0.000 


83.3 ± 2.4 ± 9.3 


18.56 ± 0.53 ± 2.1 


4.562 ± 0.040 ± 0.50 


1 


,229 ± 0.018 ± 0, 


,13 


0, 


,000 


- 0.125 


86.6 ± 1.7 ±9.7 


18.14 ± 0.53 ± 2.0 


4.535 ± 0.037 ± 0.50 


1 


,245 ± 0.017 ± 0, 


,14 


0. 


,125 


- 0.250 


84.6 ± 1.7 ± 9.5 


18.58 ± 0.60 ± 2.1 


4.526 ± 0.037 ± 0.50 


1 


.218 ± 0.017 ±0, 


,13 


0, 


,250 


- 0.375 


84.5 ± 1.7 ±9.5 


18.63 ± 0.60 ± 2.1 


4.673 ± 0.040 ± 0.51 


1, 


.262 ± 0.018 ± 0, 


,14 


0, 


,375 


- 0.500 


81.8 ± 1.8 ± 9.2 


17.53 ± 0.63 ± 1.9 


4.531 ± 0.041 ± 0.50 


1, 


.242 ± 0.019 ± 0, 


,14 


0, 


,500 


- 0.625 


78.0 ± 1.8 ± 8.8 


18.52 ± 0.58 ± 2.1 


4.262 ± 0.040 ± 0.47 


1, 


.159 ± 0.019 ±0, 


,13 


0. 


,625 


- 0.750 


70.1 ± 1.7 ± 7.9 


14.96 ± 0.51 ± 1.7 


3.820 ± 0.039 ± 0.42 


1 


.005 ± 0.017 ±0, 


,11 








5.00 - 5.50 


5.50 - 6.50 


6.50 - 8.00 




8.00 - 10.00 










(pb/(GoV/c)^) 


(pb/(GoV/c)^) 


(pb/(GeV/c)^) 




(pb/(GeV/c)=) 


-0, 
— 


,750 
,625 


0.625 

- -0.500 


175.4 ± 8.4 ± 19 
210.7 ± 7.4 ± 23 


31.5 ± 2.1 ± 3.5 
43.7 ± 2.1 ± 4.9 


2.06 ± 0.25 ±0.24 


0, 


.108 ± 0.049 ±0, 


,013 


-0, 
-0, 


,500 
,375 


0.375 

0.250 


250.9 ± 8.1 ± 28 
303.5 ± 8.9 ±34 


51.5 ± 2.3 ±5.8 

64.6 ± 2.6 ± 7.2 


4.42 ± 0.34 ±0.51 


0, 


.039 ± 0.024 ±0, 


,005 


-0, 
-0, 


,250 
,125 


0.125 

- 0.000 


323.5 ± 8.7 ±36 
352.5 ± 8.9 ± 39 


72.7 ± 2.8 ±8.1 
75.0 ± 2.7 ±8.4 


6.76 ± 0.43 ±0.79 





,197 ± 0.066 ± 0, 


,024 


0, 

0, 


,000 

,125 


- 0.125 

- 0.250 


377.2 ± 8.9 ±42 
379.9 ± 8.8 ± 42 


87.6 ± 2.8 ±9.8 
81. 7± 2.7 ± 9.2 


8.35 ± 0.45 ±0.97 





,351 ± 0.079 ± 0, 


,043 


0, 
0, 


,250 
,375 


- 0.375 

- 0.500 


383.9 ± 9.1 ±42 
371.3 ± 9.5 ± 41 


82. 4± 2.8 ± 9.2 
71.2 ± 2.8 ± 8.0 


7.33 ± 0.45 ± 0.85 





.389 ± 0.092 ±0, 


,047 


0, 
0, 


,500 
,625 


- 0.625 

- 0.750 


340.1 ± 9.3 ± 38 
296.4 ± 8.8 ± 33 


75.3 ± 2.9 ± 8.4 
56.3 ± 2.5 ± 6.3 


6.42 ± 0.48 ± 0.75 


0, 


.132 ± 0.068 ±0, 


,016 
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TABLE IV: The averaged invariant differential cross section {Eda/d^p) as a function of rapidity and pr for the inclusive 
reaction 7r~p — » tt^X at 515 GeV/c. 





Vcr 






PT (GoV/c) 












2.50 - 3.00 


3.00 - 3.50 


3.50 - 4.00 




4.00 - 4.50 








(nb/(GeV/c) = ) 


(nb/(GeV/c)=^) 


(nb/(GeV/c) = ) 




(nb/(GcV/c)2) 


-0, 
-0, 


,750 - 
,625 - 


-0.625 
-0.500 


288 ± 70 ± 33 


51 ± 15 ±5.7 


5.5 ± 3.7 ±0.6 


1 
2 


.95 ± 0.29 ± 0.21 
.92 ± 0.34 ± 0.32 


-0, 
— 0. 


,500 - 
.375 - 


-0.375 
—0.250 


297 ± 62 ± 34 


50.1 ± 9.3 ± 5.6 


3.8 ± 2.9 ± 0.4 


3 
3 


.55 ± 0.34 ± 0.39 
.79 ± 0.28 ± 0.42 


— 0. 
— 0. 


.250 - 
,125 - 


-0.125 
0.000 


291 ± 31 ± 33 


76 ± 12 ± 8.6 


16.0 ± 3.0 ± 1.8 


3 
3 


.56 ± 0.29 ±0.39 
.85 ± 0.21 ± 0.42 


0, 

0. 


,000 - 
,125 - 


0.125 
0.250 


346 ± 30 ± 40 


76 ± 10 ± 8.5 


22.5 ± 3.9 ± 2.5 


4 
4 


.59 ± 0.22 ± 0.51 
.03 ± 0.20 ± 0.44 


0. 
0. 


,250 - 
375 - 


0.375 
0.500 


367 ± 25 ± 42 


67.2 ± 9.2 ± 7.5 


20.6 ± 3.3 ± 2.3 


4 
3 


.19 ± 0.20 ± 0.46 
.96 ± 0.22 ± 0.44 


0, 
0, 


,500 - 
.625 - 


0.625 
0.750 


286 ± 27 ± 33 


72 ± 11 ± 8.1 


16.6 ± 4.5 ± 1.8 


3 
3 


.86 ± 0.25 ± 0.42 
.29 ± 0.23 ± 0.36 








4.50 - 5.00 


5.00 - 5.50 


5.50 - 6.50 




6.50 - 8.00 








(nb/(GeV/c)^) 


(pb/(GeV/c)^) 


(pb/(GeV/c)^) 




(pb/(GcV/c)2) 


-0, 


.750 - 


-0.625 


0.488 ± 0.090 ± 0.054 


118 ± 29 ± 13 


45 ± 14 ±5.0 




3.2 ± 1.9 ±0.4 


— 0. 


,625 — 


—0.500 


n «7 -1- 1 1 + 07 


259 ± 67 ± 29 


68± 18 ± 7.6 




-0, 


500 - 


-0.375 


0.89 ± 0.10 ±0.10 


316 ± 54 ±35 


75± 16 ± 8.4 




2.2 ± 1.3 ±0.3 


-0. 


,375 - 


-0.250 


0.892 ± 0.094 ±0.098 


275 ± 45 ± 30 


70 ± 17 ±7.8 




-0, 


250 - 


-0.125 


1.116 ± 0.095 ± 0.12 


362 ± 49 ± 40 


61 ± 14 ± 6.8 




7.4 ± 2.5 ±0.9 


-0. 


125 - 


0.000 


1.15 ± 0.10 ± 0.13 


404 ± 54 ± 45 


61 ± 13 ± 6.9 




0. 


,000 - 


0.125 


1.105 ± 0.094 ± 0.12 


322 ± 45 ± 36 


88 ± 16 ±9.8 




9.2 ± 2.7 ± 1.1 


0, 


,125 - 


0.250 


1.43 ± 0.11 ± 0.16 


338 ± 48 ± 37 


68 ± 13 ± 7.7 




0, 


.250 - 


0.375 


1.274 ± 0.099 ±0.14 


377 ± 50 ± 42 


102 ± 17 ± 11 




7.5 ± 2.4 ±0.9 


0, 


.375 - 


0.500 


1.093 ± 0.098 ±0.12 


357 ± 51 ± 39 


58± 15 ± 6.5 




0, 


,500 - 


0.625 


1.08 ± 0.11 ±0.12 


250 ± 47 ± 28 


56 ± 14 ±6.3 




7.6 ± 2.7 ±0.9 


0, 


,625 - 


0.750 


0.92 ± 0.10 ±0.10 


220 ± 46 ± 24 


57± 16 ± 6.4 





TABLE V: Invariant differential cross section {Eda/d^p) per 
nuclcon for the inclusive reaction 7r~Be — > r7X at 515 GeV/c, 
averaged over the rapidity interval —0.75 < j/cm < 0.75. 



PT 


TT^Bo @ 515 GoV/c 


(GeV/c) 


(nb/(GcV/c)^) 


3.00 


- 3.20 


54.5 ± 7.5 ± 6.9 


3.20 


- 3.40 


23.8 ± 3.9 ± 3.0 


3.40 


- 3.60 


16.9 ± 2.1 ± 2.1 


3.60 


- 3.80 


7.91 ± 0.97 ±0.94 


3.80 


- 4.00 


4.42 ± 0.20 ± 0.52 


4.00 


- 4.20 


2.440 ± 0.099 ± 0.28 


4.20 


- 4,40 


1,506 + 0,002 + 0.17 


ipb/(Gc\-/<.:)^) 


4.40 


- 4.60 


922 ± 41 ± 107 


4.60 


- 4.80 


571 ± 25 ± 66 


4.80 


- 5.00 


353 ± 16 ± 41 


5.00 


- 5.25 


192.9 ± 9.5 ± 22.4 


5.25 


- 5.50 


110.4 ± 6.3 ± 12.9 


5.50 


- 5.75 


69.3 ± 4.7 ±8.1 


5.75 


- 6.00 


33.4 ± 3.0 ±3.9 


6.00 


- 6.50 


20.1 ± 1.5 ± 2.4 


6.50 


- 7.00 


6.23 ± 0.73 ±0.75 


7.00 


- 8.00 


1.21 ± 0.22 ± 0.15 


8.00 


- 9.00 


0.224 ± 0.083 ± 0.029 


9.00 


- 10.00 


0.035 ± 0.039 ±0.005 



TABLE VI: Invariant differential cross section [Eda/d^p) for 
the inclusive reaction 7r~p —> r/X at 515 GeV/c, averaged over 
the rapidity interval —0.75 < j/cm < 0.75. 



PT 


TT^p @ 515 GeV/c 


(GcV/c) 


(nb/(GeV/c)^) 


3.00 - 3.50 


69 ± 62 ± 9 


3.50 - 4.00 


9.1 ± 3.8 ± 1.1 


4.00 - 4.50 


1.76 ± 0.31 ± 0.20 


(pb/(GeV/c)-) 


4.50 - 5.00 


400 ± 110 ± 50 


5.00 - 5.50 


121 ± 30± 14 


5.50 - 6.00 


19 ± 17 ±2 


6.00 - 7.00 


7.8 ± 4.7 ±0.9 


7.00 - 8.00 


0.37 ± 0.97 ±0.05 
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TABLE VII: The averaged invariant differential cross section (^Eda/dPp) per nucleon as a function of rapidity and pr for the 
inclusive reaction 7r~Be — » rfS. at 515 GeV/c. Units are pb/(GeV/c)^. 







PT (GeV/c) 






3.00 - 4.00 


4.00 - 4.50 


4.50 - 5.00 



-0.750 0-625 680± 260 ± 79 287 ± 88 ± 33 

-0.625 Q 5Q0 14000 ± 6200 ± 1700 §80 ± 220 ± 100 278 ± 61 ± 32 

-0.500 0.375 1290 ± 180 ± 150 339 ± 51 ± 39 

-0.375 - -0.250 16900 ± 5600 ± 2100 i8io±170±210 456 ±45 ±53 

-0.250 - -0.125 1710 ±150 ±200 580 ±53 ±67 

-0.125 - 0.000 27100 ± 4500 ± 3300 2210 ± 140 ±260 556±47±64 

0.000 - 0.125 „_-„n-i- Tsnn-i-c/inn ^''^0='= ^^"='=2°° 578 ± 46 ± 67 

0.125 - 0.250 2'°°" 2^°" •^^O" 2400 ± 140 ± 280 650 ± 55 ± 75 

0.250 - 0.375 „„„„„ , „ .„„ , „Q„„ 2360 ± 150 ± 270 683 ± 51 ± 79 

0.375 - 0.500 23700 ± 2400 ± 2900 228O ± 140 ± 260 688 ± 51 ± 80 

0.500 - 0.625 1 ai nn 4- 01 nn 4- oonn "^^^^ =*= =*= "^^^ 650 ± 46 ± 75 

0.625 - 0.750 2^°" 22°" 2070 ± 120 ± 240 595 ± 43 ± 69 

5.00 - 5.50 5.50 - 6.50 6.50 - 8.00 



-0.750 - -0.625 37 ±15 ±4.3 6.8 ± 4.6 ± 0.8 

-0.625 - -0.500 79 ±22 ±9.1 22.9 ±4.8 ±2.7 

-0.500 - -0.375 141 ±22 ±16 18.5 ± 4.3 ± 2.2 

-0.375 0.250 125 ± 17 ± 15 32.2 ± 5.6 ± 3.8 

-0.250 0.125 157±17± 18 36.5±4.5±4.3 

-0.125 - 0.000 143 ±19 ± 17 42.5 ±5.5 ±5.0 

0.000 - 0.125 137±19± 16 33.4±5.2±3.9 

0.125 - 0.250 197 ±22 ± 23 45.6 ± 5.5 ± 5.4 

0.250 - 0.375 225 ±21 ±26 48.0 ± 5.4 ± 5.6 . ,« , „ ar 1 n 

0.375 - 0.500 218 ±22 ± 25 46.0 ± 6.4 ± 5.4 4.36 ± 0.85 ± 0.53 



0.48 ± 0.48 ±0.06 
1.31 ± 0.48 ± 0.16 
3.63 ± 0.72 ± 0.44 
3.57 ± 0.78 ±0.43 



0.->0n n.fi2"> 181 I 20 I 21 11.-) I .^..8 I 1.9 

u.(i2.-) u.rrio i7:i±iK±2o ;i.-).9 i ri.i ± 4.2 



t 0.47 



TABLE VIII: The averaged invariant differential cross section (Eda/dPp) as a function of rapidity and px for the inclusive 
reaction 7r~p — > rjX at 515 GeV/c. Units are pb/(GeV/c)^. 



!/cT 




PT (GcV/c) 






4.00 - 5.00 


5.00 - 6.00 


-0.75 - 


-0.50 


550 ± 550 ± 64 




-0.50 - 


-0.25 


950 ± 460 ±110 


74 ± 42 ± 8.6 


-0.25 - 


0.00 


1170 ± 360 ± 140 


39 ± 38 ± 4.5 


0.00 - 


0.25 


1840 ± 430 ± 210 


147 ± 47± 17 


0.25 - 


0.50 


1020 ± 270 ± 120 


153 ± 50 ± 18 


0.50 - 


0.75 


960 ± 220 ± 110 


44± 43 ± 5.1 



